Shocked ammonia in the WR nebula NGC 2359 
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ABSTRACT 

We report the detection of the (1,1) and (2,2) metastable lines of ammonia (NH3) in the 
molecular cloud associated with the Wolf-Rayet (WR) nebula NGC 2359. Besides the CO and 
H2, this is the first molecule detected in the environs of a WR star. Width (AV1/2 = 3kms~^) 
and radial velocity (T^lsr ~ 54kms~^) indicate that the NH3 lines arises from the molecular cloud 
which is interacting with the WR star. The rotational temperature derived from the (1,1) and 
(2,2) line intensity ratios is about 30 K, significantly larger than the typical kinetic temperature 
of the ambient gas of ~ 10 K. The derived NH3 abundance is ~ 10"^. Linewidth, abundance and 
kinetic temperature can be explained if NH3 is released from dust grain mantles to the gas phase 
by shocks produced by the expansion of the bubble created by the WR stellar wind. We briefly 
discuss the implications of the detection of warm NH3 associated with a WR star in connection 
to the hot NH3 emission detected in the Galactic Center and in the nuclei of external galaxies. 
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Stars: individ- 



The interstellar medium surrounding Wolf- 
Rayet (WR) stars is expected to be strongly dis- 
turbed by the fast stellar evolution with rapidly 
varying mass loss rates, luminosities, and elemen- 
tal composition of the ejecta. The Wolf-Rayet 
(WR) star HD 56925 -WR 7 in the catalogue of 
van der Hucht (2001)- excites the wind-blown 
bubble NGC 2359 (Chu, Treffers & Kwitter 1983). 
The chemical enrichement (Esteban et al. 1990), 
the nearly spherical morphology and the kine- 
matics of NGC 2359 suggest a relatively recent 
origin for this nebula, probably produced in the 
WR stage of the exciting star. The Hii region 
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which surrounds the wind-blown bubble seems to 
be confined by molecular material in the southern 
direction (Schneps et al. 1981). The spatial dis- 
tribution of the 1-0 (1) line of H2 (St-Louis et al. 
1998) suggests that the H2 vibrationally excited 
emission is produced by the interaction of the WR 
star with the molecular material. 

Recently Rizzo, Martin-Pintado & Mangum 
(2001, paper I) have mapped the CO J = 1 ^ 
and J = 2 — > 1, and the ^'^CO J = 1 ^ emission 
over the whole nebula. They found that the two 
velocity components toward the southern part of 
the nebula have LSR radial velocities (Vlsr) of 
54 and 67 km s~^. The 67 km s~^ component 
has a narrow linewidth (1 km s^^) indicating that 
this component traces quiescent material. On the 
other hand, the component with Vlsr of 54 km 
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is broad (5 km s^^), surrounds the southern 
and eastern parts of NGC2359 (see Fig. la), and 
is spatially anticorrelated with the narrow compo- 
nent. The broad linewidths and the spatial dis- 
tribution of the 54 km s~^ component indicate 
that it is affected by the interaction of the WR 
star with the ambient material at 67 km s~^. The 
shocks produced by the expansion of the stellar 
wind bubble are likely the origin of the line broad- 
ening. 

Molecular chemistry can be an important tool 
to firmly establish if the molecular material in the 
broad component is affected by the strong shocks 
produced by the stellar wind from HD 56925. In 
particular, the ammonia (NH3) fractional abun- 
dance changes from 10"^" in photodissociation re- 
gions to 10~^ in shocked regions. This is due to the 
fact that the NH3 is a very fragile molecule, easily 
photodissociated by UV radiation (van Dishoeck 
1988; Roberge et al. 1991), but it can also be 
released from icy grain mantles by shocks with 
moderate velocities (Flower & Pineau-Des-Forcts 
1994). In this letter we report the first detection 
of ammonia emission towards a WR nebula. The 
rather high abundance of NH3 and the kinetic tem- 
perature found in the broad velocity component is 
in agreement with the idea that the molecular gas 
has been shocked by the stellar wind from the WR 
star HD 56925. 

2. Observations and results 

The (1,1) and (2,2) metastable ammonia lines 

were observed using the Effelsberg 100~m radiote- 
lescope of the MPIfR in December 2000 and March 
2001. The half power beam width of the tele- 
scope at the rest frequency of the NH3 lines, 23.7 
GHz, was 40". We used the new cooled dual chan- 
nel HEMT K-band receiver with a typical system 
temperature of 230 K on a main beam brightness 
temperature scale. An 8192-channel autocorrela- 
tor was used as the backend. Both ammonia lines 
were observed simultaneously with a total band- 
width of 20 MHz and a velocity resolution of 0.06 
km s~^. In order to improve the S/N ratio, the 
spectra were later smoothed to a resolution of 0.5 
km s~^. The observations were made in a dual 
beam switching mode (switching frequency 1 Hz) 
with a beam throw of 121" in azimuth. Point- 
ing was regularly checked and it was found to be 



within 5". The data were calibrated using the con- 
tinuum emission of NGC 7027 as the principal cal- 
ibrator, assuming a flirx density of 5.4 Jy at 23.9 
GHz (Ott et al. 1994). Orion-KL (Hermsen et al. 
1988) was used as a line calibrator. Both calibra- 
tions agree to within 10% and we conclude that 
line intensities on a main beam temperature scale 
are accurate to within ±15%. 

Figure la shows the spatial distribution of the 
"'^^CO J = 2 ^ 1 emission of the broad compo- 
nent at Vlsr = 53 km s~^ in NGC 2359, overlaid 
on an optical image of the WR nebula. The filled 
square indicates the location of the maximum CO 
emission, which is also the position where the NH3 
data were taken. Figs, lb to le show the line pro- 
files of the J = 1 ^ hues of ^^CO and "CO, and 
the (1,1) and (2,2) lines of NH3, respectively. All 
the lines were measured with similar beam sizes 
of 40"-50". Both ammonia lines are weak but 
clearly detected at a 7-8 a levels. The observa- 
tional parameters of the lines have been derived 
from gaussian fits to the profiles of the (1,1) and 
(2,2) lines, and the results are shown in Table 1. 
The radial velocities and linewidths are in good 
agreement with those of the broad CO component 
at 53-54 km s~^, indicating that both CO and 
NH3 emission arises from the same region. 



Tabk" 1: Obsc^iAatioiial and plu sical parameters 



Parameter 


Value 


RA, Dec (B1950) 


07'"16™19",-13°11'26" 


NH3 (1,1) fit: 




Peak emperature 


53 ± 6 mK 


LSR velocity 


53.1 ±0.6 km s"^ 


Width profile 


2.6 ±0.6 km s"^ 


NH3 (2,2) fit: 




Peak emperature 


46 ± 7 mK 


LSR velocity 


53.6 ±0.6 km s'^ 


Width profile 


3.0 ±0.7 km s-i 


N(CO) 


~ 7.0 X 10^6 cm-2 


N(NH3) 


~ 6.8 X 10^2 


X(NH3) 


~ 1.0 X 10"^ 


Trot 


~ 30 K 



2 



If3 
05 

S 

d 

3 -13°10'00" 
to 

g 

(1) 




7''16"40'' 20" 00' 

Right Ascension (B1950) 




0.04 



0.04 



20 40 60 80 

LSR velocity (km s 

Fig. 1.— 12 CO, "CO and NH3 emission toward 
the WR nebula NGC 2359. The intensity distribu- 
tion map in a) corresponds to the CO J = 2 — !■ 1 
emission (see Paper I). The individual profiles are 
from the black square plotted onto the clump to 
the southeast of the optical nebula. The position 
of the WN4 star HD 56925 is indicated by an ar- 
row. In the profiles from b) to e), the molecular 
transition is indicated at the top right corner. 



Assuming optically thin emission for both am- 
monia lines, we have derived the (1,1) and (2,2) 
column densities given in Table 1. Prom the ra- 
tio of the (1,1) and (2,2) column densities wc es- 
timate a rotational temperature of 30 K. Under 
typical conditions of the molecular gas, this ro- 
tational temperature can be considered lower 
limit to the kinetic temperature of the gas (Walm- 
sley & Ungerechts 1983). Using a kinetic temper- 
ature of 30 K and considering only the five lowest 
metastable (J = K) levels for the partition func- 
tion, we derive a lower limit to the total ammo- 

-2. For a CO 
(paper I) and 
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nia column density of 7 x 10^^^ 
column density of 0.7 x 10^^^ cm 
assuming a fractional CO abundance of 10~^, we 
derive a NH3 abundance of 1 ±0.4 x 10"^. This 
NH3 abundance is comparable to those measured 
in cold galactic molecular clouds (Benson & My- 
ers 1983) and shocked regions (Martfn-Pintado & 
Cernicharo 1987; TafaUa & Bachiller 1995), but 
significantly higher than in PDR's (Fuente et al. 
1990). 

3. Discussion 

The broad widths and the radial velocities of 
the (1,1) and (2,2) NH3 lines indicate that the NH3 
emission arises from the gas which is interacting 
with the WR nebula. Two different kinds of in- 
teraction between an evolved massive star and the 
surrounding material are expected: UV photodis- 
sociation and shocks produced by the stellar wind. 
It is well known that the UV radiation from hot 
stars dissociates fragile molecules, such as NII3, 
even up to visual extinctions of 5 mag (i.e., H2 
column densities of 5 x 10^^ cm"^). In this case, 
the photodissociation rate would be of the order 
of 10"'^ s-i (van Dishoeck 1988). Hence, the rel- 
atively high NII3 abundance is striking in view 
of the rather low H2 column density of ~ 10^^ 
cm~^ derived from the CO data. If the total 
extinction between the star and the NH3 region 
would be comparable to the total extinction de- 
rived from the CO data, the NII3 must be pro- 
duced very efficiently in order to compete with its 
photodissociation due to the WR star radiation 
field. Prom the dust continuum emission measured 
by the IRAS HiRcs data at 60/xm and lOO/im (Au- 
mann, Fowler & Melnyk 1990; J. R. Rizzo et al., 
in preparation), we have derived a dust temper- 
ature of 29 K and a column density in the range 
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2 — 5 X 10"'^^ cm~^ towards the position where NH3 
is detected. These values were obtained from the 
color-corrected 60/im and 100/xm emission maps 
smoothed to an angular resolution of 2', and as- 
suming an emissivity proportional to A~^. The 
dust grains were assumed to be spherical of size 
0.1/zm and with a density of 3 gcm^'^ (Hildcbrand 
1983). Assuming the standard gas-to-dust ratio of 
100, the total gas column density associated to the 
warm dust at 29 K is only 10^^ cm~^, at least two 
orders of magnitude below the total gas column 
density derived from CO in paper I. This sug- 
gests that basically all the dust associated with 
the molecular gas can not have temperatures of 
> 10 K. The warm dust component measured by 
IRAS is likely associated to the material affected 
by the UV radiation from the WR star, shielding 
the warm gas from dissociation. The gas in the 
broad component seems to be warmer than the 
dust. 

The presence of shockfronts acting in this re- 
gion could explain both the NH3 abundance and 
the gas kinetic temperature higher than the dust 
temperature in the broad component. A C-shock 
of 10-15 km would heat the gas up to the tem- 
peratures derived from our data and release the 
NH3 from the dust grain icy mantles (Flower & 
Pineau-Des-Forets 1994) to achieve the observed 
abundance. Garcia-Segura & Mac Low (1995a, 
1995b) have modelled the dynamical evolution of 
WR bubbles and predicted the presence of sev- 
eral shockfronts running into the surrounding gas 
at the different stages of the evolution of massive 
stars. When the WR stellar wind shocks the previ- 
ous RSG ejecta, an expansion velocity of ~ 10-15 
km s~^ is expected for the shocked gas. This ve- 
locity is in good agreement with the velocity sep- 
aration of 12-13 km s~^ measured between the 
shocked and the quiescent gas (paper I). 

The Galactic center region (GC) is known to 
contain large amounts of warm gas with relatively 
high abundance of molecules like NH3, SiO and 
C2H5OH extended over relatively large regions 
(Hiittemeister et al. 1995; Martfn-Pintado et al. 
1997; Martin-Pintado et al. 2001). This warm gas 
seems to be associated with the cold dust, and a 
shock chemistry has been claimed to explain the 
heating, the abundances and the densities of these 
molecules in the GC. However, the origin of the 
widespread shocks is, so far, unknown. Martin- 



Pintado et al. (1999) have recently found that a 
large fraction of the hot NH3 in the Sgr B2 enve- 
lope is located in shells of 1-2 pc, expanding at 
8-10 km s~^. The large NH3 abundance and the 
rather high kinetic temperature found in the Sgr 
B2 envelope can be explained in terms of a C-shock 
chemistry (Flower, Pineau Dcs Forcts & Walmslcy 
1995). Martfn-Pintado et al. (1999) proposed that 
these hot molecular shells are the shocked layers 
produced by the strong winds of WR stars in the 
Galactic Center. The detection of warm NH3 as- 
sociated with the expanding shells produced by a 
WR stellar wind supports the possible origin for 
the warm NH3 in the Sgr B2 envelope. Although 
the densities in the GC clouds are certainly higher 
than the density in NGC 2359, the ammonia abun- 
dances and the expansion velocities in the GC 
shells are comparable to those found in NGC 2359. 

Hot ammonia has also been observed in the nu- 
clei of the external galaxies IC 342 (Martin & Ho 
1986; Ho et al. 1990) and Maffei 2 (Henkel et al. 
2000; Takano et al. 2000). The rotational tem- 
peratures of 50 K and 80 K for IC 342 and Maffei 
2, respectively, are similar to those observed in 
the GC. In both galaxy nuclei, the derived NH3 
abundances are also similar to the GC. Further- 
more, the gas temperatures derived from NH3 are 
higher than the dust temperature and the nucleus 
of IC 342 has overall properties similar to those of 
the GC region (Downes et al. 1992). In both galac- 
tic nuclei, IC 342 and Maffei 2, the gas heating is 
not due to recently formed stars but to large scale 
shocks. As in the case of the GC, the origin of such 
shocks remains unknown. Based on high angular 
resolution data (5") from the nucleus of IC 342, 
Ho et al. (1990) conclude that the morphology of 
the NH3 emission may be a distributed bump of 
small and hot molecular clouds surrounding young 
OB stars. However, the NH3 is hard to be heated 
by the OB field radiation since it is easily pho- 
todissociated by the UV photons. It is worth not- 
ing that this ammonia distribution would also be 
consistent with hot shells like those observed by 
Martin-Pintado et al. (1999) in the GC. Further 
studies of the chemistry, physical conditions and 
kinematics of the molecular gas in the environs of 
galactic evolved massive stars would provide un- 
valuable constraints to the models for the impact 
of the evolution of massive stars onto the interstel- 
lar medium. 
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